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ABSTRACT 


A comparative study of the A hyperon equ atio ns of state of Banik , Hempel 
and Banyopadhyay (BffB) (IBanik et al. 2014) and S he n et ah ( 20111 ) (denoted 
as ffShen A) for core collapse supernova (CCSN) simulations is carried out in this 
work. The dynamical evolution of a protoneutron star (PNS) into a black hole 
is investigated in core collapse supernova simulations in the general relativistic 
one dimensional code using the BHBA</> and HShen A equation of state (EoS) ta¬ 


b les a nd different progenitor models from Woosley and Heger flWooslev & Heger 


20071 ). Radial profiles of the mass fractions of baryons, the density as well as 


the temperature in the PNS at different moments in time, are compared for 
both EoS tables. The behaviour of the central density of the PNS with time 
is demonstrated for those two A hyperon EoS tables and compared with their 
corresponding nuclear EoS tables. It is observed that the black hole formation 
time is higher in the BHBA(/> case than in the HShen A EoS for the entire set 
of progenitor models adopted here, because the repulsive A-A interaction makes 
the BHBA0 EoS stiffen Neutrino emission with the A hyperon EoS ceases earlier 
than that of its nuclear counterpart. The long duration evolution of the shock 
radius and gravitational mass of the PNS after a successful supernova explosion 
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with enhanced neutrino heating are studied with the BHBA0 EoS and s20WH07 
progenitor model. The PNS is found to remain stable for 4 s and might evolve 
into a cold neutron star. 

Subject headings: equation of state - stars:black holes - starsmeutron - super¬ 
novae: general 


1. Introduction 


The recent discovery of a 2 M 0 neutron star puts stringent conditions on the com position 


and equation of state (EoS) of dense matter in neutron star interiors flAntoniadis et ah 
2013). It has been observed that the presence of strangeness degrees of freedom such as 


hyperons makes the EoS softer, which is inco mpatible with the massive neutron star in mos t 


cases. This is known as the hyperon puzzle (IBuballa et ah 112014 iLonardoni et ah 1120151) . 
Describing hyperon matter in neutron stars is a challenge in many-body theories. It has 
been argued that the hyperon-hyperon repulsive interaction due to the exchange of strange 
vector meson makes the EoS stiffer and might overcome the puzzle. 

The /^-equilibrated equations of state that include hyperons were constructed after 


state are found to resi 

lit in 2 Mm or heavier neutron stars ( 

Weissenborn et al. 

2012a 

b; 

Lastowiccki et al. 

2012 

Coined & Sedrakian 2013 

LoDes & Menezes 

2013; 

Gusakov et ah 

2014; 

van Dalen et ah 

2014; 

Char & Banik 

2014) 

Besides hyperons, the antikaon conden- 


sate w as also included in some calculations, which led to massive neutron stars flChar fe Banik 
20141 ). In all of these calculations, the repulsive hyperon-hyperon interaction that is mediated 
by (/> mesons was considered. 

Many EoS tables involving hypero ns were dev eloped for supernova simulations. The 
first hyperon EoS table was prepared by iTshiznka et al. (120081). In this case, the full baryo n 


octet was added to the Shen nucleon EoS table fllshizuka et al. l2008t IShen et ah 19981 ). 


The Shen nucleon EoS table was based on a relativistic mean held (RMF) model that had 
the Thomas-Fermi appro ximation for the description of inhomogeneous matter below the 
saturation density flShen et al. 1119981) . Later, the Shen nucleon E oS was extended to include 
only A hyperons in the HShen A EoS table flsiien et ah 2011 ). Another extensively used 
supernova EoS is the Lattimer and Swesty (LS) nucleon EoS table, which based on the 
non-relativistic Skyrrne interactio n flLattimer fc Swestv 19911 ). Recently A hyperons were 
included in the LS nucleon EoS (lOertel et al. 1120121) . All these hyp eron EoS tables were 


used in core collapse supernova (CCSN) simulations by several authors fllshizuka et al. 2008 
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Nakazato et al. II20081.120121; ISnmivoshi et al. Il2009t iBanik II20131.12014) . However, none of 


these hyperon EoS tables were consistent with the 2 M e neutron star constraint. 

Recently, we computed EoS tables that included A hyperons within the fr ame work 
of the density dependent relativistic hadron (DDRH) held theory flBanik et al. 112014) . In 
those EoS tables, light and heavy nuclei, as well as interacting nucleons are described in 
the n uclear statistical equilibrium model which takes into a ccount the excluded volume ef¬ 
fects flHempel fe Schaffner-Bielich 1 12010 c Banik et al. 2014). Two variants of the hyperon 
EoS tables were generated - in one case (BHBA0), the repulsive A hyperon - A hyperon 
interaction mediated by (f) mesons was considered, and in the other case (BHBA), this inter¬ 
action was neglected. It should be noted that the DDRH model with the DD2 parameter 
set for nucleons is in very good agr e emen t with the symmetry en ergy properties at the sat¬ 


uration density flTvpel et al. 2010. 2013; Lattimer & Lim 120131) . We imposed the charge 


neutrality and ^-equilibrium conditions on the BHB hyperon EoS tables and calculated the 
mass-radius relationship of the neutron star sequence. It was observed that the maximum 
mass corresponding to the BIID Ao EoS w as 2.1 M 0 which is well above the recently observed 
massive neutron star flBanik et al. 2014). Other hyperon EoSs for /3-equilibrated neutron 
star matter gave r ise to the maximum mass n eutron stars of 1.75 M m for the HShen A EoS 
(Slicn et. al. 2011 ) , 1.6 M m for Ishizuka EoS f lshizuka et al. 2008 ). 1.91 M 0 for the LS+A 
EoS f Oertel et al. 2012 ). 


I 11 this paper, for the first time, we perform supernova simulations with the BHBA</> 
EoS table, which is compatible with a 2 Mq neutron star, in the general relativistic one 
dimensional (GRID) hydrodynamics code flO’Connor fe Ott 20101 ). Our main goal is to in¬ 
vestigate the appearance of A hyperons in the postbounce phase and the role of A hyperons 
in CCSNs. Furthermore, we compare our simulation results with those of previous c alcula- 
tions with other hyperon EoS tables, particularly the HShen A EoS table flBanik 20 131). We 
are looking for important effects of hyperons in CCSN with the BHBA0 EoS compared with 
those of other hyperon EoS. 


The paper is organised as follows. In Section 2, the DDRH model for BHBA</> EoS 
table and the RMF model for the HShen A EoS table are described. We also briefly discuss 
the GRID model for CCSN simulations. The results of our calculation are discussed and 
compared with those of the HShen A EoS in Section 3. Section 4 contains the summary and 
conclusions. 
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2. Methodology 


First First we shal 


discuss the salient feature of the BH BAc/> and HShen A EoS tables 


for CCSN simulations flShen et al. Il201ll; IBanik et ah I2014T) . The EoS tables are functions 
of three parameters i.e. baryon number density, temperature, and proton fraction. In both 
cases, the compositions of matter that vary from one region to the other depending on 
those parameters are nuclei, (anti)neutrons, (anti)protons, (anti)A hyperons, photons plus 
electrons and positrons that form a uniform background. The contribution of (anti)neutrinos 
is not added to the EoS tables and is dealt with separately. We describe the baryonic 
contribution below. 


2.1. BHBA0 and HShen A EoS tables 


In the BHBA0 EoS table, the nuclear statistical equilibrium (NSE) model of lHem pel fe Schaffner -Bielic 


(1201 Oil is applied for the description of the matter made of light and heavy nuclei, and un¬ 
bound nucleons at low temperatures and below the saturation density, whereas the high 
density matter is described within the framework of the DDRH model adopting the RMF 
approximation (IBanik fe Bandvopa,dhvavll2002l; ITvpel et al. II2010I; IBanik et ah 11201411 . The 
repulsive interaction between A hyperons mediated by 0 mesons is included in the RMF 
model. Nucleon-meson couplings in the DDRH model are density dependent. The DD2 
parameter set of nucleon-meson couplings i s used to describe the nucl ear matter properties 


(ITypel fe Wolter Ill999l: ITvpel et ah 120101.120131: iFischer et al. 120141) . It should be noted 


t hat the n uclear EoS in the DDRH model using DD2 parameter set is known as HS(DD2) 
( Fischer et al. 2014 1. 


On the other hand, the uniform matter at high density and temperature in the HShen A 
EoS table was descri bed wit hin the framework of the RMF model including nonlinear terms in 
cr and oj mesons (IShen et ah 201111 ; non-uniform matter at low temperatures and below the 
saturation density was considered as a mixture of alpha particles, heavy nuclei, and unbound 
nucleons. Heavy nuclei were calculated using the Thomas-F ermi approac h. The_Shen EoS 
exploited the single nucleus approximation for heavy nuclei (IShe n et al. 19981. 201111. The 
interaction among A hyperons due to </> mesons was neglected in this case. Furthermore, in 
this case baryon-meson couplings of the RMF model are density-independent. We denote the 
EoSs wit h and w ithout A hype rons as HShen A and HShen, respectively. The parameter set 
from Sugahara & Tok i (1199411 that is known as the TM1 set was adopted for the nucleon- 
meson coupling constants of the RMF model. 


The nuclear matter saturation properties of two RMF models discussed above are 
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recorded in Table Q] It should be noted that though the incompressibility of nuclear matter, 
symmetry energy, and its slope coefficient of t he DD2 set at the saturation density are i n 
very good agreement with experimental values (ILat tim er fe L im 2013; Fischer et ah 2014), 
the corresponding quantities of the TM1 set are not. This would h ave ser ious bear ing on 
the description of high density matter in the RMF model of HShen (IShen et al. 2 01111 . For 
both EoS tables, A hyperon-vector meson couplings are estim ated from the SU(6) symmetry 


relations flDover fe Gal 11985c ISchaffner fe Mishustin 1 1199611 and A hyperon - scalar meson 


coupling is obtained from the hypernuclei data. The A hyperon potential depth is -30 Me V 
in normal nuclear matter (IMillener et ah 1 11988c iMares et ah Il995l; ISchaffner et ah 119921) . 


The EoSs of /3-equilibrated and charge neutral cold neutron star matter with and with¬ 
out A hyperons are calculated from the supernova EoS Tables. The maximum masses of 
cold neutron stars without A hyperons for HS(DD2) and HShen EoS are given by Table [U 
Furthermore, the maximum masses of cold neutron stars corresponding to the BHBA0 and 
HShen A are 2.1 M 0 and 1.75 M 0 ( Banik et ah 2014 ; Shen et ah 2011 ). respectively. 


For CCSN simulations, we make use of the HS(DD2), BHBAA HShen and HShen A 
EoS tables which are available from the stellarcollapse.org website □. 


2.2. General relativistic model for supernova simulations 

We perform the CCSN simulations using the sphericall y symmetr ic gene ral relativistic 
hydrodynamics code GRID which was developed by 10’Connor fe Ott (120101) . Microphysi¬ 
cal EoSs for supernova matter and an approximate treatment of neutrinos in the pre- and 
postbounce phases are implemented in the GRID code. We use the BHBA0 and HShen 
A EoS tables in CCSN simulations with the GRID code. Three neutrino species d enoted 
by is e , u e , v x {— Vfx, v T , v T ) are considered in this model (jO’C onn or & Ott 120111) . Key 
aspects of neutrino heating and cooling are incorporated into the model. The leakage 


scheme (IRuffert et al. Il996l; iRosswog fe Liebendorfer 120031) exploited in the GRID code 


gives approximate number and energy emission rates. The neutrino heating rate consid¬ 
ered here involves the scale factor fheat which could be enhanced beyond the normal value 


of 1 to achi eve ad ditional neutrino heating for ’’successful” CCSN explosions (Jankaj 2001 


O’Connor fe Ott 2 0111 ). We take fheat = 1 hi CCSN simulations, if not stated otherwise. 


In principle, an accurate and expensive neutrino treatment should be based on the 
Boltzmann neutrino transport. However, computationally efficient schemes for neutrinos 


^ee http://stellarcollapse.org/equationofstate 
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are employed in the GRID code for faster CCSN simulations. Moreover, it lias been noted 
that the results obtained in CCSN simulations using the simplified treatment of neutrino 
leakage and heating in the GRID were quantitatively similar to the results obtained from 
one dimensional (ID) simulations with the Boltzmann neutrino transport by other groups 
(iFischer et al. 1120091; ISumivoshi et al. Ir2009l i. It was argued that progenitor structures played 
more important roles in the coll apse o f a protoneutro n star (PNS) to a black hole than the 
details of neutrino treatment flO’Connor fe Ott 20111). 


3. Results and Discussion 


Now we report our investigations on CCSNs within the GRID code us ing the HShen A 


hyper on and BHBA0 EoS tables. In these studies, nonrotat ing progenitors oflWooslev fe Heger 


( 120071) (WH07) are used. In their stellar evolution studies IWooslev fe Heger (120071) evolved 
zero age main-sequence (ZAMS) stars with solar metallicity denoted by the prefix s before 
presupernova models, follo wed by ZAMS mass. Significant mass loss was reported in sWH07 
presupernova models (jO’Connqr Sz Ott 20111). 


We perform the CCSN simulations with presupernova models as recorded in Table [2] 
In all numerical calculations, we fix the neutrino heating factor jheat = 1. In the next 
paragraphs, we discuss the results of simulations starting from the gravitational collapse 
of the iron core followed by the core bounce to the postbounce evolution of the PNS for 
s40WH07 and s23WH07 models with the HShen A and BHBA0 EoS tables in details. In all 
of these simulations, a shock wave is launched at the core bounce, it stalls after traversing 
a few 100 km, then recedes and becomes an accretion shock. Because neutrinos in the 
ID CCSN model could not revive the shock, the PNS shrinks due to mass accretion and 
its density and temperature increase during the postbounce evolution. This leads to the 
appearance of A hyperons in the PNS. 


For s40WH07, the core bounce occurs at 0.273 and 0.321 s, corresponding to the HShen 
A hyperon and BHBA0 EoS, respectively. Similarly, in the s23WH07 model the core bounce 
times for the HShen A and the BHBA0 EoS are 0.266 and 0.315 s, respectively. The ap¬ 
pearance of strangeness or A hyperons in the postbounce phase and its role in the evolution 
of the PNS are the main focuses of this investigation. For s40WH07 and s23WH07 models 
and both hyperon EoS tables, A hyperons do not populate the PNS at the core bounce. In 
simulations with both presupernova models, strangeness in the form of A hyperons sets in a 
few hundred milliseconds (ms) after the core bounce and increases with time thereafter. 

Figure 1 depicts the PNS compositions as a function of radius at two different postbounce 
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times for s40WH07 with the HShen A (left panel) and BHBA</> (right panel) EoS tables. For 
postbonnce time (tpb) 0.31 s, the central value of A fraction is higher for the BHBA0 EoS 
than that of the HShen A hyperon EoS. The profile of A hyperons is wider in the latter case. 
We find similar trends for A hyperons at a later time t p b = 0.51 s. For both EoS tables, 
the population of As increases with time, ft is to be noted that the central value of the A 
fraction is a high density effect, whereas the off-centre As are populated thermally. We study 
the density and temperature profiles to understand this behaviour. 

The density prohles as a function of radius are plotted for s40HW07 at the bounce as 
well as for t p b = 0.31 and 0.51 s in Figure 2. The left panel of the figure corresponds to the 
HShen A EoS and the right panel implies the results of the BHBA</> EoS. At the bounce, the 
central density (p c ) of the PNS in both cases is just above the normal nuclear matter density, 
as evident by the figure. Though the density prohles for both EoS tables are quantitatively 
the same at t p b = 0, they differ at later times. The central density at t p b = 0.51 in the right 
panel is higher than that of the left panel. In both cases, the central density exceeds two 
times the normal nuclear matter density. This high central density facilitates a significant 
population of As in the core of the PNS, as seen in Fig. 1. However, the density falls well 
below normal nuclear matter density at the tail of the profile. The off-centre As in Fig. 1 
could not be attributed to the density effect. 

The temperature prohles as a function of radius are shown for s40WH07, with the 
HShen A hyperon (left panel) and the BIIBAd (right panel) EoS tables in Figure 3. Just 
as in Figure 2, the temperatures prohles are plotted at the core bounce and t pb = 0.31 and 
0.51 s in both panels of Figure 3. The peaks of temperature prohles located away from the 
centre of the PNS for both EoSs after the core bounce later shift toward the centre with 
time in both panels. It is to be noted that the central temperature at the bounce is higher 
for the BHBA</> EoS compared with the corresponding temperature for the HShen A EoS. 
Furthermore, the peak temperature around 8 km at 0.51 s after core bounce in the case of 
the BHBA0 EOS is much higher than the corresponding scenario for the HShen A EoS. This 
high temperature results in thermally produced A hyperons away from the centre of the PNS 
as shown in Figure 1. We find from Figure 1 that thermal As are more abundant around 8 
km at later times for the BHBA0 EoS due to a higher peak temperature. 

We also study prohles of particle fraction, density, and temperature for s23WH07 using 
both hyperon EoS tables as shown in Figures 4-6. We obtain qualitatively similar results for 
s23WH07 as we have already discussed for s40WH07. 

Now we focus on the postbounce evolution of the PNS for different presupernova models 
with nuclear and A hyperon EoS tables corresponding to the HShen and BHB models. Figure 
7 exhibits the evolution of the central density of the PNS with the postbounce times for 


s40WH07 (left panel) and s23WH07 (right panel). Results are shown in both panels for the 
HShen nuclear EoS, the HShen A EoS, the HS(DD2) nuclear EoS and the BHBA</> EoS. It 
should be noted that the core bounce time for the hyperon EoS is the same as that of the 
corresponding nuclear EoS. In all cases in both panels of Figure 7, we find that the central 
density increases gradually to several times the normal nuclear matter density. Finally, 
there is a steep rise in the central density when the PNS dynamically collapses into a black 
hole in milliseconds. It should be noted that the black hole formation time is different for 
different EoS models. It is evident from the CCSN simulation of s23WH07 that the black 
holes are formed at 1.511 and 1.623 s after the core bounce for the HS(DD2) and the HShen 
EoS, respectively. For s40WH07, the black hole formation time is 0.942 s in the case of 
the HS(DD2) EoS, whereas it is 1.084 s for the HShen EoS. For both supernova models and 
nuclear EoS tables, the black hole is formed earlier in case of the HS(DD2) than the situation 
with the HShen EoS. The maximum gravitational (baryonic) PNS masses are 2.464 (2.616) 
M 0 and 2.459 (2.587) M 0 for s40SW07 with the HS(DD2) and the HShen EoS, respectively. 
Similarly, for s23WH07, those are 2.428 (2.605) M 0 and 2.383 (2.512) M 0 corresponding to 
the HS(DD2) and the HShen EoS. On the other hand, the dynamical collapse to a black 
hole is accelerated for the HShen A and BHBA0 EoS tables because hyperons make the EoS 
softer. It is evident from Figure 7 that the black hole formation time is shorter for hyperon 
EoS than that for the corresponding nuclear EoS. However, there is little difference between 
the black hole formation times corresponding to the HShen A and BHBA0 EoSs. 


The results of CCSN simulations with other presupernova mo dels are recorded in Table 
[2 The first column of the table lists the presupernova models of Woos lev & Heg er 020071 ) 
starting from s20WH07 to s40WH07. Two EoS tables, such as the HShen A and the BHBA0 
are adopted in these calculations. Under each EoS, the first column represents the black 
hole formation time ( tsH ) estimated from the core bounce and the next column gives the 


maximum baryon mass (Mb }Tnax ) followed by the maximum gravitational mass (M g}Tnax ) of 
the PNS at the point of instability corresponding to the central value of the lapse function 
0.3. Further investigations with the two A hyperon EoSs reveal an opposite behaviour of tsH 
than what has been observed for nuclear EoSs. For A hyperon EoS, tsH for the BHBA0 is 
always greater than that of the HShen A for all presupernova models except s40WH07. The 
comparison of two hyperon EoSs shows that the BHBA^ is a stiffer EoS than the HShen 
A. The explanation of this behaviour may be traced back to the inclusion of repulsive A-A 
interaction in the BHBA0 EoS. For all presupernova models and EoSs adopted in simulations, 
it is evident from the table that the maximum gravitational mass of the PNSs in each case 
is higher than their corresponding maximum cold neutron star masses. However, in some 
cases, the maximum gravitational mass of the PNS collapsing into a black hole with the 
HShen A EoS is less than the two solar mass limit because the HShen A EoS does not result 
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in a 2 M 0 cold neutron star. It is interesting to note that in the case of the HShen A EoS, 
the difference between M. g max of the PNS and the maximum mass of the cold neutron star 
that includes A hyperons (1.75 M 0 ) is appreciable, whereas the maximum gravitational mass 
of the PNS for the BHBA0 EoS is very similar to the value of the corresponding maximum 
mass of the cold neutron star with As (2.1 M 0 ) for the entire set of progenitor models. This 
shows that the thermal effects in the PNS for the BHBHA0 might not be as strong as in the 
PNS with the HShen A because the EoS is stiffer in the former case. The role_of decreasin g 


therm al pressure with increasing stiffness of the EoS was already noted by lO’Connor fe Ott 
12 0111 ). This should have interesting implications for the study of the metastability of the 
PNS with the BHBA</> EoS. 

We compare our findings with other CCSN simulations with hyperon EoS. The Ishizuka 
hyperon EoS includ es A, E, and 5 hyperons and is an extension of the HShen nuclear 
EoS (llsh i zuka e t al. 120081 ). The CCSN simulations were performed in a spherically sym¬ 
metric genera l re lativistic n eutrino radiation hydrodynamics model usin g a 40 progen- 
itor of 


Wooslev fc Weaver (119951 ) and the Is hizuka hyperon Eo S flSumivoshi et ah 


Nakaza.to et al. 2 012 ). With the LS+A EoS f Oertel et al. 11201 2) Peres et al. f Peres et ah 


2009 


20131 ) carried out a s imilar investigat i on us ing an s40W W progenitor a nd a low metallicity 


40 M 0 progenitor of IWooslev et ah 1120021) called «40. iBanik 1120131) also studied C CSN 


simulations using the HShen A EoS and progenitor models of Woosl ev & Heger 1)20071 ). par¬ 
ticularly studying the long duration evolution of the PNS in the context of understanding the 
fate of the compact object in SN1987A. It should be noted that though our results with the 
BHBA0 EoS are qualitatively similar to those of earlier calculations, they are quantitatively 
different because only our A hyperon EoS is compatible with the 2 M 0 limit of cold neutron 
stars. The early black hole formation due to softening in the A hyperon EoS compared with 
the nuclear EoS is a robust conclusion in all of these calculations. Total neutrino luminosity 
as well as v e , P e , and v x luminosities as a function of postbounce time are plotted in Figure 8 
for the HS(DD2) (left panel) and the BHBA0 (right panel) EoS. The results are shown here 
for the s40WH07 model. It should be noted that the neutrino emission ceases earlier for the 
BIIBAc) case than for the scenario with the HS(DD2) nuclear EoS. A similar c onclusion wa s 
arrived at in the simulation with other hyperon EoSs (Sumiyoshi et ah 120091 ; Banik 120131 ). 
The shorter neutrino burst corresponding to the A hyperon EoS before the collapse of the 
PNS into the black hole could be an important probe for the appearance of A hyperons in 
the PNS. This demands a more accurate treatment of neutrinos in the GRID code. Figure 9 
exhibits the neutrino luminosities for both the A hyperon EoS and the s40WH07 model. We 
find similar features for neutrino luminosities for both cases. Though we are considering a 
phase transition from nuclear to A hyperon matter, we do not find any evidence for a second 
neutrino burst, which was observed in a first order quark-hadron phase and was responsible 
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for a successful supernova explosion (ISagert et al. 1120091) . 


So far we have seen that simulations in the ID CCSN model might lead to accretion 
driven black holes in failed supernovae. If a successful supernova occurs, can exotic matter 
such as hyperons make the PNS metastable and drive it to become a low mass black hole 
during the long duration evolution when thermal support decreases and deleptonization takes 
place in the P NS? Such a scenario was envisaged for the non-ob 


in SN1987A flBrown fe Bethe 1 19041: IPrakash et al. 


This prob lem was also studied in CCSN simulations flKeil fe Janka 


& Banc 

voDadhvav 

2001) 

1995: 

Baumearte et al. 


1996 ; Banik 20131) . We continue our study by increasing the neutrino heating scale factor 


to fheat — 1-5 for s20WH07 with the BHBA0 EoS. The left panel of Figure 10 exhibits 
the shock radius as a function of postbounce time. For the neutrino scale factor fheat — 
1, it fails to launch a successful supernova explosion and the shock radius recedes. Finally, 
the PNS collapses into a black hole. For fheat = 1-5, it is observed that the shock radius 
increases with time after a successful supernova explosion. The PNS remains stable until 
4 s. We do not find any onset of the metastability in the PNS due to the loss of thermal 
support and neutrino pressure during the cooling phase over a few seconds. The window for 
the metastability is very narrow because the maximum PNS mass in this case is 2.138 M 0 
whereas the maximum cold neutron star mass corresponding to the BHBA0 EoS is 2.1 M 0 . 
The PNS might evolve into a cold neutron star. Gravitational masses of the PNS for fheat — 
1 and 1.5 are shown as a function of postbounce time in the right panel of Figure 10. The 
PNS cools down to a neutron star with a mass ~ 1.64 M 0 at the end of 4 s. 


4. Summary and Conclusions 


We have performed CCSN simulations using the BHBA0 EoS, which is co mpatib l e wit h a 
2 Mr.-, neutron star, and several progenitor models from the stellar studies of Wooslev _& Heg er 
(120071 ). It is observed that As are produced a few hundred milliseconds after the core bounce. 
The appearance of A hyperons is studied in great detail. It is evident from the density and 
temperature profiles as a function of radius that As are produced in the core of the PNS when 
the central density exceeds two times the normal nuclear matter density during the post¬ 
bounce evolution phase. On the other hand, an off-centre population of thermal A hyperons 
is the result of peak values of temperature away from the centre of the PNS. When we set 
the neutrino heating scale factor fheat — 1, each CCSN simulation ends with the formation 
of a black hole driven by mass accretion. It is interesting to find out that the black hole 
formation time for the BHBA0 EoS is shorter than that of the HShen A EoS though the 
opposite conclusion is drawn from the accretion driven black hole with the HShen nuclear 
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and HS(DD2) EoS models. This is attributed to the fact that the repulsive A-A interaction 
in the BHBA0 EoS makes it a stiffer EoS than the HShen A EoS. Neutrino luminosity is 
found to cease with the formation of a black hole earlier for the A hyperon EoS than for 
the corresponding case with the nuclear EoS. We have studied the nretastability of the PNS 
due to the BHBA0 EoS in the long duration evolution after a successful supernova explosion 
using the s20WH07 progenitor model with the increased neutrino heating scale factor of 
fheat = 1-5. In this case, we do not find any delayed collapse into the black hole due to the 
presence of A hyperons in the PNS. The PNS that has a mass ~ 1.64 M e remains stable 
until 4 s and might become a cold neutron star. 


The numerical calculations presented in this article have been performed in the blade 
server of the Astroparticle Physics and Cosmological Division, Saha Institute of Nuclear 
Physics. We acknowledge fruitful discussions of GRID with Evan O’Connor. DB thanks the 
Alexander von Humboldt Foundation for support. 
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Table 1. The saturation properties of nuclear matter such as saturation density (no), 
binding energy (BE), incompressibility (K), symmetry energy (S), and slope coefficient of 
symmetry energy (L) are obtained using the DD2 and TM1 parameter are obtained using 

the DD2 and TM1 parameter sets. 


Parameter 

n 0 

BE 

K 

S 

L 

M max 

Set 

( fmr 3 ) 

(MeV) 

(MeV) 

(MeV) 

(MeV) 

(M s ) 

DD2 

0.1491 

16.02 

243 

31.67 

55.04 

2.42 

TM1 

0.1455 

16.31 

281 

36.95 

110.99 

2.18 


Note. — Maximum masses of cold neutron stars without A hy- 
perons corresponding to the HS(DD2) and the HShen EoS are also 
mentioned here. 
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Table 2. Black hole formation time, baryoni c and gravitatio nal masses of PNSs for CCSN 
simulations with the progenitor models of Wooslev fe Heger (120071 1 and the BHBA0 and 

HShen A EoS tables. 


Model 

t-BH 

(*) 

BHBA0 

Mb,max 

(Me) 

Mq max 

(M e ) 

t-BH 

(s) 

HShen A 

Mb,max 

(Me) 

Mg,max 

(Me) 

S20WH07 

1.938 

2.251 

2.138 

1.652 

1.999 

1.964 

S23WH07 

0.879 

2.276 

2.203 

0.847 

2.095 

2.073 

S25WH07 

1.548 

2.234 

2.141 

1.376 

2.035 

2.001 

S30WH07 

2.942 

2.243 

2.113 

2.258 

1.967 

1.929 

S35WH07 

1.175 

2.243 

2.161 

1.084 

2.071 

2.041 

S40WH07 

0.555 

2.250 

2.210 

0.565 

2.129 

2.118 


Note. — For all cases considered here, fheat = 1- 
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Fig. 1.— Mass fractions of different species in the PNS are shown as a function of radius 
for the HShen A EoS (left panel) and the BHBA0 EoS (right panel) at £ p & = 0.31 and 0.51 
s (online-version: red). The results in both panels correspond to the s40WH07 model. 
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Fig. 2.— Density profiles of the PNS are shown as a function of radius for the HShen A EoS 
(left panel) and the BHBA0 EoS (right panel) at the core bounce and t p b = 0.31 (online- 
version: red) and 0.51 s (online-version: green). The results in both panels correspond to 
the s40WH07 model. 
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Fig. 3.— Temperature profiles of the PNS are shown as a function of radius for the HShen A 
EoS (left panel) and the BHBA</> EoS (right panel) at the core bounce and t pb = 0.31 (online- 
version: red) and 0.51 s (online-version: green). The results in both panels correspond to 
the s40WH07 model. 
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Fig. 4.— Same as Figure 1 but for the s23WH07 model. The results correspond to the HShen 
A EoS (left panel) and the BHBA</> EoS (right panel) at t p b = 0.31 and 0.51 s (online-version: 
red). 
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Fig. 5.— Same as Figure 2 but for the s23WH07 model. The results correspond to the 
HShen A EoS (left panel) and the BffBA0 EoS (right panel) at the core bounce and t p & = 
0.31 (online-version: red) and 0.51 s (online-version: green). 
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Fig. 6.— Same as Figure 3 but for the s23WH07 model. The results correspond to the 
HShen A EoS (left panel) and the BHBA0 EoS (right panel) at the core bounce and t p b = 
0.31 (online-version: red) and 0.51 s (online-version: green). 

























22 




Fig. 7.— Central baryon density is plotted with the postbonnce time for the HShen nuclear 
EoS, the HShen A EoS (online-version: red), the HS(DD2) (online-version: green) and the 
BHBA<^> EoS (online-version: blue). The results in the left and right panels correspond to 
the s40WH07 and s23WH07 models. 
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Fig. 8.— Total neutrino luminosity (online-version: blue) as well as v e , u e (online-version: 
red) and u x (online-version: green) luminosities are plotted with the postbounce time for 
the HS(DD2) (left panel) and the BHBA0 (right panel) EoS. The results correspond to the 
s40WH07 model. 
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Fig. 9.— Total neutrino luminosity (online-version: blue) as well as u e , v e (online-version: 
red) and v x (online-version: green) luminosities are plotted with the postbounce time for 
the HShen A (left panel) and the BHBA0 (right panel) EoS. The results correspond to the 
s40WH07 model. 
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Fig. 10.— Shock radius (left panel) and gravitational mass of the PNS (right panel) are 
plotted with the postbounce time using the neutrino heating factor f heat = 1 and 1.5 (online- 
version: red) for the s20WH07 model and the BHBA0 EoS. 














